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ABSTRACT 


Illicit drug-trafficking is a major concern of the United States and is a primary pillar of President 
Barack Obama’s Strategy to Combat Transnational Organized Crime. In the eastern Pacifie and 
Caribbean Sea, drug-trafficking organizations operate a variety of vessels to transit drugs from 
South America to the United States. Joint Interagency Task Force (JIATF) South, in cooperation 
with partner agencies and nations, detects, tracks, and interdicts illegal drug-trafficking in this 
region. In this thesis, we develop a probability model based on intelligence inputs to generate a 
spatial temporal heat map specifying the likely location of targets over time. We also formulate a 
path-finding model that takes the heat map as input and determines route characteristics through 
transit, such as departure times, waypoints, and speed. We link the results of our models to a 
separate effort that seeks to provide JIATF South with an optimal search plan to maximize the 
expeeted amount of drugs seized. We show that our path-finding model aceurately reproduces 
the target’s track even when the target transits along a eomplex route. Furthermore, we show 
that the optimal search plan based on our path-finding model is nearly identical to the search 
plan using the known parameters. 
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Executive Summary 


Illicit drug-trafficking originating from South America is a major concern of the United States, 
and continued expansion of drug-trafficking significantly affects the national security interests 
of the United States. Organizations tasked with combatting drug traffickers are faced with 
several problems. They must detect, track, and interdict illegal drug traffickers en route to 
various follow-on transit zones in Central America, Mexico, and the Caribbean. They utilize 
incoming intelligence from various sources to determine where and what assets to employ. 
Their efforts, while productive, are not able to stop a large majority of the drugs from reaching 
the United States. They battle Mexican and Colombian Drug Trafficking Organizations (DTOs) 
spread throughout Latin America. The DTOs have extensive capabilities and access to large 
amounts of money. Combined with the close proximity to the United States, the DTOs pose an 
extreme threat to the nation’s security. 

Operators at Joint Interagency Task Force (JIATF) South obtain various types of intelligence 
information and utilize that information when creating search and interdiction plans. They may 
have at their disposal information concerning departure times, departure and arrival locations, 
type of vessel, speed, drug load, and waypoints. At all times, there is a level of uncertainty as¬ 
sociated with the intelligence, which can be represented by a probability distribution for each of 
the aforementioned pieces of information. Given this information and surrounding uncertainty, 
JIATF South operators can benefit from a visual heat map displaying the likely route of a target 
or targets over time. In this thesis, we address the problem of creating a heat map of likely drug 
trafficker locations based on intelligence inputs. 

One approach to generating a heat map is through simulation. Simulations have the benefit of 
being able to account for many real world complexities like weather. However, they also may 
be computationally expensive and difficult to validate and verify. We desire to complement 
and vet such complicated simulation by developing more analytically tractable models. We 
develop probability models to generate a heat map to provide operators with better situational 
awareness about the likely locations of targets, both in the current period and future periods. 
We first build a simple probability model that accounts for uncertainty only in the departure 
and arrival locations. Using this model as a baseline, we further enhance its functionality to 
account for uncertain speed and time of departure, as well as consideration for multiple legs, 
multiple targets in the same region, and the ability to output expected amount of drugs. Given 
the probability distributions for the intelligence inputs (e.g., speeds, locations, and times), the 




model produces the spatial temporal probability map of the trafficker’s path. Using the heat 
map as input, we then address the problem of determining a drug trafficker’s transit route based 
on the collected intelligence information. 

After developing a model to transform intelligence inputs into a heat map, we next turn to the 
inverse problem: Given a heat map, what are the most likely intelligence inputs used to create it. 
We adopt this approach to compare the path-finding model’s results with some known inputs via 
the search vehicle optimization model. Additionally, it is possible that the only data available 
for JIATF South planners will be a heat map, thus providing a means to estimate intelligence 
parameters is beneficial for the planners. 

We produce the necessary intelligence parameters through the use of various regression tech¬ 
niques and consideration for the input requirements of the optimization model. The path-finding 
model approximates the path taken by the target as a sequence of piecewise linear legs, where 
the collection of all legs produces the target’s entire track from departure location to arrival 
location. Standard regression techniques are incapable of automatically identifying where a 
new leg begins, which occurs when speed or heading changes. As a result, we employ more 
advanced regression techniques, namely Multivariate Adaptive Regression Splines (MARS). In 
accordance with the search vehicle optimization model, we define a leg as a segment of the 
target’s track that has a constant heading and speed. This creates a piecewise linear path for the 
target. Specifically, we determine the target’s departure time range, expected departure time, 
arrival time range, expected arrival time, speed, departure and arrival locations, and waypoint 
locations. 

The final aspect of the path-finding model is specifying the Area of Uncertainty (AOU) of a tar¬ 
get as it traverses along a path leg. Roughly speaking, the AOU is the size of the target’s "blob" 
on the heat map. The AOU along a target’s track represents the uncertainty in both departure 
speed and departure time. With these uncertainties, an area exists in which it is possible for the 
target to be present at any point in the associated area. The search plan algorithm requires that 
the AOU be constant on each leg of the path. Thus, if the AOU changes significantly, we break 
one leg of the path into multiple legs even if the speed and heading do not change. We develop 
some heuristics to accomplish this. 

We develop two case studies and the corresponding intelligence information from which we 
generate a heat map according to the probability model. Next, we use our path-finding algo¬ 
rithm to estimate the input parameters and compare how well they match to the true values. 


xviii 



We show the eapability to identify speed ehanges, ehanges in heading or eourse, and ehanges 
in the AOU surrounding the target throughout the path. We illustrate the models’ funetionality 
with two ease studies of varying eomplexity. First, a simple straight path between the origin 
and destination shows the path-finding model’s ability to identify the target’s starting loeation, 
destination, speed, and assoeiated travel times. Combining these route eharaeteristies, we are 
able to reproduee the target’s likely traek from start to finish. Supplying our path-finding results 
into the seareh vehiele alloeation model shows similar seareh plan results when eompared with 
the original ease parameters used in the probability model. Next, we demonstrate a more eom- 
plex seenario with a target traveling to the destination with a eourse ehange during the journey. 
Our model aeeurately identifies this eourse ehange and reproduees the target’s traek. These two 
ease studies are representative of actual drug trafficker routes but additional analysis is needed 
to validate the model’s ability. Finally, we provide the search vehicle allocation model with the 
path-finding model results for a scenario with two targets. The resulting search plan is nearly 
identical to the one generated when the true intelligence inputs are known. This provides strong 
evidence that we have provided a key link between the probability, path-finding, and optimiza¬ 
tion models. 
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CHAPTER 1: 

Introduction 


Illicit drug-trafficking is a major concern of the United States and is a primary pillar of Pres¬ 
ident Baraek Obama’s Strategy to Combat Transnational Organized Crime. Transnational Or¬ 
ganized Crime (TOC) eneompasses many areas of global erime to inelude human smuggling, 
traffleking in persons, weapons traffieking, intelleetual property theft, eybererime, and drug- 
traffieking (National Seeurity Staff, 2011). Drug-traffieking has expanded in recent years, and 
organizations have beeome inereasingly more violent and less eentralized. The U.S. govern¬ 
ment lists Drug Traffieking Organizations (DTOs) as the "greatest organized erime threat" to 
the United States (National Drug Intelligenee Couneil [NDIC], 2008). Murder rates in known 
drug-traffieking eountries have inereased, and eontinued expansion of drug-traffieking signif- 
ieantly affeets the national seeurity interests of the United States. This thesis eontributes to 
the ongoing researeh to improve the effeetiveness of the ageneies tasked with eombatting il- 
lieit drug-traffieking. Based on ineoming intelligenee, the analysis models movements of drug 
traffiekers. 


1.1 Background 

Latin Ameriea is a prominent player in the global drug market, with involvement in nearly 
all areas of the illieit drug supply ehain. Aetivities in the region inelude drug erop farming, 
drug production, drug-trafficking, and drug consumption. The South Ameriean eountries in the 
Andean region, partieularly Colombia, Peru, and Bolivia, are the sole produeers of eoeaine for 
the global market (United States Department of State, 2012). United States Department of State 
(2012) estimates total eoeoa produetion in these three eountries at 187,500 heetares in 2010. 
This is down from a peak in 2007 of 232,500 heetares. Combining the Andean region’s ability 
to produee eoeaine and their elose proximity to the United States, makes Latin Ameriea the 
eentral hub of the illieit drug market (Seelke, Wyler, Beittel, & Sullivan, 2011). Additional 
faetors that eontribute to the region’s role in the drug market inelude poverty, inequality, the 
absenee of other fiseal opportunities, and a weak governing foree (Seelke et ah, 2011). These 
faetors have left DTOs nearly free to aet as they wish in the supply of illieit drugs. 

A variety of DTOs operate in Latin Ameriea, and they range from small family-based operations 
to global networks that are well eonneeted throughout the world to inelude ties with terrorist 
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organizations. The more advaneed DTOs employ eonsiderable paramilitary and eounterintelli- 
genee assets to eontend with legitimate poliee forees and eonduet operations with minimal risk 
via eorrupt governments (Seelke et ah, 2011). The family-based operations typically specialize 
in a few areas of the drug supply chain. 

Mexican and Colombian DTOs are the most prominent in the region. They currently have 
operations worldwide to include Europe, West Africa, and the Asia-Pacific regions as well as 
connections to organized crime and terrorists (Seelke et ah, 2011). The NDIC (2008) reports 
that Mexican and Colombaian DTOs generate, remove, and launder between $18 billion and $39 
billion in wholesale drug proceeds annually. Currently, Mexican DTOs are the dominant groups 
supplying most illicit drugs to the United States. Seven Mexican DTOs are currently operating 
in all nine of the Organized Crime Drug Enforcement Task Eorces (OCDETE) regions. Table 1.1 
and Eigure 1.1 show the Mexican DTOs and their primary drugs and regions of operations. 



Eigure 1.1: Nine regions of the Organized Crime Drug Enforcement Task Eorce. All seven 
Mexican-based DTOs operate in all of the OCDETE regions. (Erom NDIC, 2011) 


Over the last ten years, the Mexican DTOs have replaced many of the Colombian DTOs that 
controlled the market in the 1980s and 1990s. The Mexican organizations have a geographical 
advantage over the Colombian drug traffickers. They have access and control of the smuggling 
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DTO 

Primary Drugs 

Primary Regions 

Sinaloa Cartel 

Cocaine 

Heroin 

Marijuana 

MDMA 

Methamphetamine 

Florida/Caribbean 
Great Fakes 
Mid-Atlantic 

New England 

New York/New Jersey 
Pacific 
Southeast 
Southwest 

West Central 

Los Zetas 

Cocain 

Marijuana 

Florida/Caribbean 
Great Fakes 
Southeast 
Southwest 

Gulf Cartel 

Cocaine 

Marijuana 

Florida/Caribbean 

Mid-Atlantic 

New England 

New York/New Jersey 
Southeast 
Southwest 

Juarez Cartel 

Cocaine 

Marijuana 

Great Fakes 

New York/New Jersey 
Pacific 
Southeast 
Southwest 

West Central 

BLO 

Cocaine 

Heroin 

Marijuana 

Southeast 

Southwest 

LFM 

Cocaine 

Heroin 

Marijuana 

Methamphetamine 

Southeast 

Southwest 

Tijuana Cartel 

Cocaine 

Heroin 

Marijuana 

Methamphetamine 

Great Fakes 

Pacific 

Southwest 


Table 1.1: Mexican-based DTOs operate in all nine of the OCDETF regions and supply most 
illicit drugs to the United States. (From NDIC, 2011) 
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routes along the U.S. southwest border, and these routes or plazas are fiercely contested among 
the Mexican DTOs (NDIC, 2011). Additionally, the ability to obtain, transport, and distribute 
the illicit drugs in demand in the United States allows them to dominate the market. The Sinaloa 
Cartel, with operations in all nine OCDETF regions, is the principal Mexican DTO. It boasts 
an extensive transportation and distribution network, a vast network of personnel to aid in traf¬ 
ficking operations, and control of heroin, marijuana, and methamphetamine production (NDIC, 
2011). According to NDIC (2011), the Mexican DTOs will remain the dominant U.S. illicit 
drug suppliers for the near future. 

The Colombian DTOs, the Cali cartel and Medellin cartel that dominated the trade in the 1980s 
and 1990s, have been replaced by Norte del Valle DTO. Rather than being one central cartel 
such as their predecessors, the Norte del Valle consists of many smaller cells and is less hi¬ 
erarchical and more decentralized (Seelke et ah, 2011). State Department designated Foreign 
Terrorist Organizations (FTOs), the Revolutionary Armed Forces of Colombia (FARC), the 
National Fiberation Army (EFN), and the United Defense Forces of Colombia (AUC), control 
and finance Colombian illicit drug-trafficking operations. The FTOs then use the drug money 
to attack Colombian police and military (United States Department of State, 2012). These or¬ 
ganizations are becoming more involved in drug-trafficking operations. They tax cocoa farmers 
and buyers and have controlling interests in various operations of the drug supply chain (Seelke 
et ah, 2011). When the AUC demobilized, smaller paramilitary groups formed and increased 
their involvement in the drug market, proving to be a major challenge for law enforcement. 
Their increased presence has caused a rise in violence in both Colombian cities and rural vil¬ 
lages (United States Department of State, 2012). 

Many of the drugs, specifically cocaine, that reach U.S. soil transit from South America through 
the Central America-Mexico transit routes or via the Caribbean transit routes en route to the 
United States. The Consolidated Counterdrug Database (CCDB) documented 1,494, 1,228, and 
804 metric tons of cocaine movement through the Western Hemisphere Transit Zone in fiscal 
years 2008, 2009, and 2010, respectively (Office of National Drug Control Policy, 2013). The 
primary transportation route is through the Central America-Mexico corridor to include terri¬ 
torial waters in both the eastern Pacific and the Caribbean Sea. This drug flow accounts for 
95% of the cocaine going to the United States. More than half of these drugs transit through 
Central America first (Seelke et ah, 2011). This is a shift from the operations of the 1980s 
and early 1990s, when drugs were mainly smuggled through the Caribbean to Florida (Seelke 
et ah, 2011). In 2010, 5% of the cocaine transiting to the United States was via the Caribbean, 
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whereas in 2001, 26% was transiting via the Caribbean (United States Senate Caueus on Inter¬ 
national Nareoties Control, 2012). In September 2011, President Obama identified 22 eountries 
as major drug transit or major illieit drug produeing countries. Of the 22 identified, 19 of the 
countries are in Latin America or the Caribbean (United States Department of State, 2012). In 
the Caribbean region specifically. The Bahamas, the Dominican Republic, Haiti, and Jamaica 
are major countries used in the illegal trafficking of cocaine (United States Senate Caucus on In¬ 
ternational Narcotics Control, 2012). Approximately one half of the cocaine that flows through 
the Caribbean transits through the Dominican Republic (United States Senate Caucus on In¬ 
ternational Narcotics Control, 2012). Puerto Rico presents another major transit hub. Since 
the country is within the U.S. customs zone, regulations to get into the United States are less 
stringent, thus 70% to 80% of the cocaine moving through Puerto Rico is bound for the United 
States while the remaining cocaine is consumed within Puerto Rico (United States Senate Cau¬ 
cus on International Narcotics Control, 2012). The Caribbean presents three main corridors in 
Figure 1.2 for the illicit transport of drugs. 



Figure 1.2: Illicit Caribbean drug routes are divided into three main corridors: Central, Eastern, 
and ABC. (From United States Senate Caucus on International Narcotics Control, 2012) 

The central corridor is a primary avenue for moving cocaine from South America through Ja¬ 
maica, Haiti, the Dominican Republic, or The Bahamas into the United States. The eastern 
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corridor moves heroin and eoeaine through eastern Caribbean nations, Trinidad, and Tobago 
to Puerto Rieo and ultimately into the United States. The ABC eorridor of Aruba, Bonaire, 
and Curaeao is prevalent among eoeaine traffiekers as well (United States Senate Caucus on 
International Nareoties Control, 2012). 

In 2011, 80% of drugs destined for the United States moved via the waterways and the remain¬ 
ing 20% via aireraft (United States Senate Caueus on International Nareoties Control, 2012). 
Maritime and air transit routes from Joint Interageney Task Foree (JIATF) South are shown 
in Figures 1.3 and 1.4. DTOs using maritime Central Ameriea-Mexico transit routes utili z e 
both Caribbean and eastern Paeifie routes as evident in Figure 1.3. The air routes leaving the 
Colombia-Venezuela border avoid Colombian air spaee in an effort to evade surveillanee by 
the U.S. Government and Colombian Government as part of the Air Bridge Denial Agreement. 
Colombian offleials eite that the surveillanee and interdietion efforts foree smugglers to remain 
clear of Colombian air space (GAO, 2009). 



Figure 1.3: Maritime drug transit routes in 2011. 80% of the drugs destined for the United 
States transit via waterways. (From United States Senate Caueus on International Nareoties 
Control, 2012) 


The primary means of transport is via maritime vessels, whieh inelude go-fast boats, pleasure 
eraft, fishing vessels, and Self-propelled Semi-submersible (SPSS). The Caribbean has seen a 
significant increase in the number of go-fast boats used to transport illicit drugs. These vessels 
travel at high speeds, are small, display nearly no radar eross-seetion, and are hard to see 
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Figure 1.4: Air drug transit routes in 2011. (From United States Senate Caucus on International 
Narcotics Control, 2012) 


during the day, all of which makes them difficult to find, track, and interdict. They can carry 
as much as two metric tons of cocaine and their transits usually last one or two days (Office of 
National Drug Control Policy, 2013). In the eastern Pacific, fishing vessels typically depart from 
Colombian and Ecuadorian ports carrying several tons of cocaine to ports in Central America 
or Mexico. These operations can last up to six weeks and inject substantial amounts of cocaine 
into the Central America-Mexico transit route. Drug traffickers have recently been employing, 
to a greater extent, more subtle maritime transportation in the form of SPSSs. 

The use of SPSSs, such as the one in Figure 1.5, to transport cocaine is on the rise in the 
Caribbean and eastern Pacific. Between 2001 and 2007, 23 events involving SPSSs were doc¬ 
umented, compared to 62 SPSS events in the first three quarters of fiscal year 2008 (USCG, 
2009). SPSSs are built mainly from wood and fiberglass in the jungles of Colombia. They 
are capable of carrying up to 10 metric tons of cocaine with a range of 2,500 nautical miles. 
They range in length from 25 to 80 feet and travel up to 13 knots. A crew of 4 to 5 operates 
the vessel, which is designed to scuttle upon detection by law enforcement (USCG, 2009). Re¬ 
cently, DTOs have been improving their capabilities with the construction of fully-submersible 
vessels. In July of 2010, the Ecuador police seized a drug submarine while under construction 
in the jungles of Ecuador that is 100 feet long and capable of diving to 65 feet (Kraul, 2010). 
The switch to fully submersible submarines poses a significantly harder threat to combat. 
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Figure 1.5: SPSS as seen by U.S. Coast Guard prior to being boarded. (From USCG, 2009) 


Antidrug efforts have been a mainstay of U.S. drug poliey in Latin Ameriea and the Caribbean 
for the last 30 years, with the United States spending $14.6 billion to fight illegal drugs. The 
strategy is to stop it at the souree by destroying drug produeing erops, interrupting shipments, 
and improving soeio-eeonomie faetors in eountries prone to DTOs. The programs are funded by 
the U.S. Department of State, U.S. Ageney for International Development, and the Department 
of Defense (DoD). The DoD works with foreign eountries to train, equip, and improve eounter 
drug efforts. National Seeurity Direetive 18, signed by President Bush on August 21, 1989, and 
the Andean Initiative put eounter nareoties on the DoD’s list of eore priorities. It is now the lead 
federal ageney responsible for deteetion and monitoring of illegal drugs from foreign eountries 
into the United States. The National Defense Authorization Aets of 1991 and 1998 as amended 
through FY2011 allows the DoD to provide support to eertain eountries, mainly Colombia. 
More reeently, the Andean Counterdrug Initiative (ACI) proposed by President Clinton in 2000 
and passed through Congress inereased funding for eounter drug operations in South Ameriea. 
The primary eountries of the ACI inelude Colombia, Peru, and Bolivia (Seelke et ah, 2011). In 
fiseal year 2008, the Merida Initiative inereased funding for eounter drug operations in Mex- 
ieo and sinee then efforts have shifted from Colombia to Mexieo, Central Ameriea, and the 
Caribbean. Two reeent initiatives, the Central Ameriea Regional Seeurity Initiative (CARSI) 
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and the Caribbean Basin Security Initiative (CBSI) have brought more focus to the region and 
DoD operations specifically support these initiatives (Seelke et ah, 2011). 

1.2 U.S. Southern Command 

DoD counternarcotics policy is established by the Office of the Deputy Assistant Secretary 
of Defense for Countemarcotics and Global Threats. Two combatant commands in the west¬ 
ern hemisphere execute the mission of providing counter drug support, U.S. Southern Com¬ 
mand (SOUTHCOM) and U.S. Northern Command (NORTHCOM). NORTHCOM’s Area of 
Responsibility (AOR) covers Mexico and SOUTHCOM’s AOR covers all of Latin America 
south of Mexico. SOUTHCOM’s efforts to counter illegal drug-trafficking directly supports 
the first objective of their Command Strategy 2020, defend the United States (Fraser, 2010). 
SOUTHCOM is composed of five component commands, three joint task forces, and one direct 
reporting unit. SOUTHCOM collaborates with partner nations to combat TOC in their AOR. 
Illicit drug-trafficking is their primary focus as this is a means for TOCs to strengthen their 
influence in the region. SOUTHCOM provides support to partner nations enabling them to pro¬ 
vide better security in the region. This support is directed to interrupt the flow of illegal drugs, 
remove TOC networks, and provide assistance to remove criminal influences in lawless areas. 
As counter narcotics are a core priority of the DoD, SOUTHCOM is the lead agency in the 
AOR in carrying out this mission. Their efforts are in line with the objectives of the Strategy 
to Combat Transnational Organized Crime, the CARSI, and the CBSI (United States Southern 
Command, 2013). 

1.3 Joint Interagency Task Force South 

JIATF South is a task force of SOUTHCOM that directly performs SOUTHCOM’s mission to 
combat illicit drug-trafficking. Located in Key West, FL, JIATF South detects, tracks, and inter¬ 
dicts illegal drug-trafficking activity in the eastern Pacific, Gulf of Mexico, and the Caribbean. 
They work with partner nations, military organizations, and U.S. law enforcement agencies 
in carrying out their mission. U.S. military and partner nations’ maritime and air assets are 
employed in the detection, monitoring, and interdiction of illicit drug-trafficking in the region 
(United States Southern Command, 2013). The assets available to JIATF South can vary over 
time and have recently been in decline. Typically, U.S. Navy, U.S. Coast Guard, and partner 
nation maritime vessels patrol the Caribbean Sea, Gulf of Mexico, and the eastern Pacific. A 
U.S. Coast Guard Law Enforcement Detachment embarks on a U.S. vessel to take control dur¬ 
ing boarding operations to seize illicit drugs and apprehend suspects (United States Southern 
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Command, 2013). In the air, U.S. military, interagency, and partner nation aircraft execute de¬ 
tection and monitoring missions and direct interdiction assets to suspicious surface vessels for 
further investigation. 

In 2011, JIATF South interdicted 119 metric tons of cocaine, arrested 355 drug traffickers, 
and seized 70 aircraft and vessels, which resulted in a $7.1 billion loss of revenue for DTOs 
(Posture Statement of General Douglas M. Fraser, 2012). They projected 775 - 930 metric 
tons of cocaine to move towards the United States in 2012 and in January, JIATF South be¬ 
gan Operation Martillo with the goal of disrupting drug-trafficking along Central American 
coastal waters (Posture Statement of General Douglas M. Fraser, 2012). U.S. Naval Forces 
Southern Command, in Mayport, FL, supported Operation Martillo with seven frigates, one re¬ 
plenishment ship, and four fixed-wing aviation squadrons (Posture Statement of General John 
F. Kelly, 2013). In the first year of Operation Martillo, JIATF South and partner nations have 
disrupted 152 metric tons of cocaine, confiscated $7.2 million in bulk cash, and seized 4 aircraft, 
8 SPSSs, and 89 vessels (Posture Statement of General John F. Kelly, 2013). This success is 
closely dependent on support from partner nations who played a role in 67% of the disruptions 
associated with Operation Martillo (Posture Statement of General John F. Kelly, 2013). JIATF 
South has played a major role in the planing and execution of this operation. Based on the 
amount interdicted in 2011, approximately 85% of the cocaine is not stopped. They oversee a 
large operational area and effective and efficient employment of available assets is paramount 
to their success. 

1.4 Mathematical Modeling Approach 

This thesis addresses some of the uncertainty involved with detection and monitoring of illicit 
drug-trafficking in the Caribbean and eastern Pacific. JIATF South is the primary agency that 
is most interested in this work. Given the number of assets they have and the vast amount 
of area under their watch, their success is heavily dependent on employing available assets as 
effectively as possible. The main purpose of this overall endeavor is to provide better situational 
awareness to operators, analysts, and targeteers about the likely location of targets in the present 
and future. To meet this goal, we formulate two models. 

JIATF South faces two primary problems in this context. First, the locations of the targets 
are unknown. A probability model creates a heat map from various intelligence inputs, which 
include departure location, arrival location, and departure time. This provides a better picture 
of the overall AOR. Next, applying a path-finding algorithm provides JIATF South with a 
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probable track of the target, which will aid in the asset allocation process. Employment of their 
assets represents the second problem. After developing an estimate for the likely location of 
a target over time, where to send their assets is the next complexity to overcome. Pietz and 
Royset (2013) develop an optimization model that takes as input the output from the path¬ 
finding algorithm and creates a search plan for JIATF South. 

Based on incoming intelligence, the probability model creates a time-layered map indicating 
the probability associated with a drug trafficker being in a specific area at a specific time. As 
additional intelligence is received, the map is updated with the new information. Following the 
creation of the probability map, an algorithm determines the target’s course, speed, waypoints 
(if applicable), time of departure, and time of arrival. The optimization model uses the algorithm 
output to maximize the effectiveness of JIATF South assets with respect to the potential amount 
of drugs seized. With the connection of our probability model and path-finding algorithm to the 
optimization model, JIATF South receives an all-inclusive tactical decision aid to aid in their 
efforts to counter illicit drug-trafficking. 

1.5 Thesis Structure 

This thesis develops a method to determine a drug trafficker’s course and speed through time 
given some initial intelligence reports and links additional work being performed in support 
of JIATF South’s mission. Chapter 2 details a literature review to support the efforts of this 
study. Chapter 3 describes the probability model, key uncertainties, and results of the probabil¬ 
ity model. Chapter 4 details the target path-finding model and key inputs to the search vehicle 
optimization model. Chapter 5 presents the conclusions and recommendations of this thesis. 
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CHAPTER 2: 
Literature Review 


This thesis covers a range of topics that have been extensively researched in the past and lever¬ 
ages techniques from probability, regression, and optimization. Topics such as search theory 
can be traced back to World War II and the work of Koopman and colleagues on the Operations 
Evaluation Group (OEG). More recently, researchers at the Naval Research Eaboratory (NRE) 
have developed a probability model of likely pirate locations based on both intelligence and 
historical weather data. A major contribution of this thesis links our path-finding model to an 
HATE South asset optimization model. This research and other applicable work is reviewed to 
give context to this thesis. 

2.1 Search Theory 

The discipline of search theory originates from the need to locate enemy submarines in World 
War II. Koopman (1946) presents initial findings in this area of research focusing on various 
types of optimal search and detection techniques with the constraint to maximize use of limited 
resources. Washburn (2002) provides many mathematical techniques to determine the proba¬ 
bility of detecting a target. Rozen (2009) develops a stochastic dynamic programming model to 
optimize the operational policy of an interdiction force performing a maritime interdiction mis¬ 
sion. These works consider cases of varying intelligence about a target, where best to search, 
and when to perform interdiction, while this thesis develops a grid of probabilities overlaid on 
the AOR and then determines the target’s track. Our focus develops tools to generate the inputs 
for search models such as Pietz and Royset (2013). 

2.2 Piracy Models 

Hansen et al. (2011) develop a probability model to aid in the location of pirates. Their model 
dynamically combines real-time intelligence and weather information to provide a visual pic¬ 
ture of likely pirate locations. Intelligence such as homeports, vessel information, performance 
characteristics, pirate activity, pirate operations, and observations of likely pirate vessels are 
inputs in the model. Additionally, wind, wave, and current data is supplied creating a "dynam¬ 
ically coupled" model (Hansen et al., 2011). The model produces a prediction of the pirate’s 
track over time with the constraints of intelligence and weather. After thousands of runs are 
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performed, the model produces a probability map or "heat map" of pirate locations. This map 
is updated as additional intelligence is received. 

Building on the Hansen et al. (2011) piracy model, An et al. (2012) proposes an extension to the 
model that dynamically allocates search and interdiction assets to minimize the probability of 
a successful pirate attack. Their research focuses on allocating multiple search and interdiction 
assets against multiple targets or pirates. An et al. (2012) divide the asset allocation problem 
into two phases. First, they allocate the interdiction assets. Second, they divide the area not 
covered by interdiction assets into rectangular areas for search asset allocation. 

A pirate’s use of a small vessel is significantly impacted by winds and sea state leading to the 
development of this model. Similarly, drug traffickers operate small vessels during their transit 
from South America to the major drug transit countries making them susceptible to the same 
types of weather. The Hansen et al. (2011) model entails a high degree of complexity through 
the simultaneous analysis of weather, intelligence, and pirate behavior. This model has been 
extended to the drug trafficker problem but retains its complexity. Our goal is to develop a 
much simpler model. We analyze straight line paths between an origin and destination without 
the addition of weather information. Expanding this approach we apply similar techniques to 
target paths that involve multiple courses and non-constant velocity. Our approach focuses on 
determining a single target’s path through the AOR. 


2.3 Probability Model - Optimization Model Link 

Pfeiff (2009), Gift (2010), and Bessman (2010) develop attacker - defender type interdiction 
models to maximize the searcher’s probability of detecting an illicit drug trafficker. These 
theses make use of intelligent smugglers who are able to learn and adapt their voyage to the 
searcher’s actions. This thesis, while retaining some elements of search and detection, is pri¬ 
marily concerned with determining the path of the target to include start time, origin, speed, 
course, waypoints (if applicable), end time, and destination. 

Pietz and Royset (2013) and Pietz (2013) develop an optimization model that routes JIATF 
South detection and interdiction assets to maximize the amount of drugs interdicted. This model 
does not take a probability map directly as input, but takes piecewise linear paths to describe 
the target track and an Area of Uncertainty (AOU) surrounding the path. JIATF South planners 
may not have access to the necessary parameters to run the optimization model and if weather 
significantly impacts the target’s course, the inputs will be insufficient. Our focus on transform- 
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ing the probability map into a path links the probabilities generated from gathered intelligenee 
to the optimization model of JIATF South assets. Furthermore, we develop an algorithm to 
determine the optimal AOU surrounding the target. These efforts directly support the inputs 
needed for the optimization model. 

2.4 Statistics Problem 

The probability model produces a set of probabilities indexed over time and identified by the lat¬ 
itude and longitude of the respective grid cell center point. To determine the target’s track during 
the time period, several statistical methods are considered. Principle curves are smooth curves 
that pass through the middle of the data cloud minimizing the distance from the data (Hastie, 
1984). These smooth curves give a general idea of the target path but lack information needed 
on specific waypoints. Smoothing splines provide similar descriptions of the target track as in 
principle curves but also lack the needed granularity on when a target changes course (Hastie, 
Tibshirani, & Friedman, 2009). Hastie et al. (2009) discusses the complexity of determining 
the appropriate location of knots, and suggests the use of Multivariate Adaptive Regression 
Splines (MARS). 

MARS employs a greedy algorithm and various approximations to reduce the complexity of 
knot selection (Hastie et al., 2009). MARS has been employed in areas of research con¬ 
cerning survival analysis, time series, ecology, and economics since its introduction by Fried¬ 
man (1991). MARS is a non-parametric regression technique that automatically models non- 
linearities and selects variables based on their contribution to overall model fit. This results in 
a piecewise linear model (Hastie et al., 2009). For this thesis, we desire to define a target’s 
track as a set of piecewise linear paths which is what MARS produces automatically. MARS 
allows paths with multiple courses or speeds to be modeled accurately by automatically de¬ 
termining locations of potential course or speed changes. Automatic selection of these points 
makes this technique preferred over other regression techniques. Once MARS determines the 
points of course or speed changes, all points are checked for validity and any spurious points are 
removed. MARS supplies the points of speed and heading change and the algorithm builds an 
ordinary least squares regression model utilizing the selected points of course or speed change. 

2.5 Area of Uncertainty 

When working with probabilities, there is an inherent uncertainty involved. In our case, this 
uncertainty resides in the location of the target. At each time step, a set of probabilities are 
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created over a particular area. An AOU is calculated for each time period using several methods, 
which include minimum bounding box, convex hull, and a greedy algorithm. 

Minimum bounding boxes and convex hulls are used in a variety of applications, including col¬ 
lision avoidance, hidden object determination, and shape analysis. Johnson and Cohen (1998), 
Lempitsky, Kohli, Rother, and Sharp (2009), and Huebner, Ruthotto, and Kragic (2008) discuss 
these applications in detail. For our purposes, we present a two-dimensional set of points that 
we wish to enclose with the minimum bounding box. The area of the minimum bounding box 
then represents the maximum AOU of the target at a particular time and encloses the entire 
probability mass. The convex hull further reduces the area while still capturing all of the prob¬ 
ability mass. When the maximum AOU or entire probability mass is not necessary, we propose 
the use of a greedy algorithm to determine the AOU for a user specified level of probability. 
This approach allows the user to select the minimum amount of probability the AOU should 
contain. 
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CHAPTER 3: 
Probability Model 


JIATF South processes various amounts of intelligenee information and makes deeisions on 
where and when to employ available assets. Given this intel, we develop probability models 
to generate a heat map to provide the operators with better situational awareness about the 
likely loeations of targets, both in the eurrent period and future periods. We first build a simple 
probability model that aeeounts for uneertainty only in the departure and arrival loeations. From 
this model, we further enhanee its funetionality to aeeount for uneertain speed and time of 
departure as well as eonsideration for multiple legs. We first present the simplest model and 
then deseribe the embellished model. 


3.1 Fundamental Probability Model 

The fundamental one target analytie model aeeounts for uneertainty only in the departure and 
arrival loeations on a reetangular area of operation. This setup does not eonsider the geography 
of realistie shorelines and we assume the target travels from the bottom of the reetangle to the 
top of the reetangle in a straight line as shown in Figure 3.1. We further assume the target’s 
departure time and veloeity are known. We set up the AOR by defining an area from 0 to 
Xmax units in width and 0 to Ymax units in height. Within the X^ax by Ymax area, we diseretize 
the AOR into a set of square grid eells for future probability ealeulations. We further assume 
that the target’s y eoordinate at the origin is 0 and the y coordinate at the destination is Tmax- 
Figure 3.1 illustrates the model setup. 

The random variables in this ease apply to the departure and arrival loeations in the x diree- 
tion and are defined as D and A respeetively. The target veloeity and departure time remain 
eonstant and are deterministie in this simplified model. Additionally, given the departure and 
arrival loeations, we assume the target moves in a straight line between the loeations. Let v 
represent the target veloeity and to represent the target departure time. The random variable D 
follows a probability distribution Fo on the interval [D/,D„] and the random variable A follows 
a probability distribution Fa on the interval [A/, A^]. Let a represent a realized value of A and d 
represent a realized value of D. We desire to find the probability that the target is in a speeifie 
box at a speeifie time. We define the box by upper and lower eoordinates in both the x and y 
direetions. Let xi and Xu define the lower and upper edges of the box in the x direetion and yi 
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Simple Probability Model Layout 



X 

Figure 3.1: Initial layout of the fundamental probability model, which shows the various depar¬ 
ture and arrival locations, target track, and box of interest. The target travels from the lower left 
comer of the area of interest to the upper right comer with constant velocity. 


and yu define the lower and upper edges of the box in the y direction. Let {xt,yt) be the coordi¬ 
nates of the target at time t. The target position, (xf,y/), is dependent on a, d, v, t, and to and is 
defined in Equations 3.1 and 3.2. We establish an indicator variable equal to 1 if the target is in 
the box and 0 otherwise. Equation 3.3 defines the indicator variable. Eor simplicity, we define 
square grid cells or boxes for probability calculations but any arbitrary shaped area can be used 
by modifying Equation 3.3. 
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if xi<xt< Xu and yi<yt< yu 
otherwise 


(3.3) 


Given a departure loeation, d in [Di^Du] and an arrival loeation, a in [A/,Ak], we know preeisely 
whether the target is in the box at a speeifie time from Equation 3.3. The probability based on 
an uneertain d and a results from integrating over the distributions, Fd and Fa, whieh yields 
Equation 3.4. 


r^u r^u 

P[xi<Xt<Xu,yi<yt<yu]= j / I{d,a,v,t,tQ)fD{d)fA{a)AdAa (3.4) 

J Di JAi 

3.2 Numerical Illustrations 

After deriving a method to determine the probability a target is in one partieular box at a speeifie 
time, we expand this to eaeh time period and all of the boxes eomprising the operating area. 
Speeifieally, for each time period, we want to compute the probability a target is in any grid cell 
making up the operating area. We present a numerical example to illustrate the application of 
the model over a multiple hour time period covering the entire area of operations. 

The simple one target analytic model as implemented in R, uses discrete departure and arrival 
distributions. We divide the departure location into rid bins and the arrival location into ria bins. 
Eor example, the uniform departure distribution takes the value 

di = T); + [Du — Di )——with probability — for 1 < z < rid (3.5) 

nd-l fid 

and the arrival random variable. A, takes on similar quantities defined as aj. Thus, the integrals 
in Equation 3.4 become sums as shown in Equation 3.6. 
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nd na 

P[xi <xt< Xu,yi <yt< yu] = IE I{di,aj,vj,xi,yi,Xu,yu)fDidi)fAiaj) (3.6) 

i=lj=l 

For example, let D and A be random variables with uniform probability distributions on the 
intervals and [A;,A„] respectively. Then Equation 3.6 becomes 

P[xi <Xt <Xu,yi<yt <yu] = - Y^Y^I{di,aj,vd,xuyi,XuAu) (3.7) 

nd na 

The probability model produces target tracks from each di to each aj and plots the associated 
tracks. When the indicator variables are summed and multiplied by the associated density func¬ 
tions, we produce the probability that the target is in a specific grid cell for each time period. By 
repeating the calculations for each time period over the entire length of the journey, we produce 
a matrix of probabilities. This matrix of probabilities is then used to generate a spatial temporal 
probability map for each time period. An example of the heat map is shown in Figure 3.2 where 
the departure location follows a uniform distribution, on the interval [10,30] and the arrival 
location follows a uniform distribution. Fa, on the interval [70,90]. 


3.3 Enhanced Probability Model 

The addition of uncertainty in velocity and departure time further enhances the fundamental 
probability model. The additional uncertainty introduces two more random variables. We de¬ 
fine these as V for velocity and T for departure time. The random variable V follows a prob¬ 
ability distribution Fy on the interval [V],K] and the random variable T follows a probability 
distribution Fr on the interval [T/, r„]. Equation 3.4 gains two additional integrals to account for 
the additional uncertainty as shown in Equation 3.8. 


Plxi <xt< Xu,yi <yt< yu] = 

rTu rVu rDu rAu 

/ / / / I{di,aj,v,t,xi,yi,Xu,yu)fT{t)fv{v)fD{d)fA{a)AtdvAdAa (3.8) 

Jt, Jv, Jd, Ja, 

Furthermore, the enhanced probability model can be used to generate probabilities for targets 
with multiple legs. The legs are divided by waypoints or geographic positions where the target 
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Figure 3.2: Heat maps from several periods of the one target simple analytie model utilizing 
uniform departure and arrival distributions on the intervals [10,30] and [70,90], respeetively. 

ehanges eourse, speed, or both. This funetionality of the enhaneed model is an extension of 
the previous model where rather than having just two points on the path, departure and arrival. 
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there can can be an arbitrary number of waypoints. If there are M total waypoints in addition to 
the departure and arrival points, then there will be M additional integrals. 

As in the fundamental probability model, the basic unit considered by the enhanced model 
is probability and if the mass of the target is all on the water, the sum of the mass is one. The 
enhanced model still utilizes square grid cells for probability calculations but the area of interest 
considers the actual geography of the eastern Pacific and Caribbean. Figure 3.3 shows the 
enhanced model’s area of interest. By considering realistic geographies Equation 3.5 becomes 
more complex in that the current 1-dimensional equations become 2-dimensional equations. 



Figure 3.3: Depiction of the area of interest in the enhanced probability model. The enhanced 
probability model considers the realistic geography of the eastern Pacific and Caribbean. 


We construct a model that generates target location probabilities over the course of the target’s 
voyage. The model takes as inputs: a starting location, a destination, target speed, waypoints, 
expected drug load, and Intel certainty. The model applies a predetermined probability distribu¬ 
tion to the locations, speed, and departure time inputs. The user can choose either a determinis¬ 
tic distribution, uniform distribution, or triangular distribution. The additional considerations in 
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the enhanced model introduce more complexity than the basic model in Section 3.1 but the core 
framework of the model is the same. Numerically, implementation of the enhanced model only 
requires additional for loops and retains its simplicity when compared to the simulation models 
such as the piracy model discussed in Section 2.2 

The enhanced model provides several generalizations to produce a more realistic heat map. The 
uncertainty in velocity and departure time simply add two integrals to Equation 3.4 and the 
more complex area of interest requires a more careful bookkeeping of departure and arrival 
locations. By including the expected drug load as an input, the model can produce a heat map 
in terms of pounds of drugs vice probability by simply multiplying the various probabilities 
by the expected drug load. Considerations for target existence can be made by multiplying the 
probabilities by the intel certainty. We can model one target or several targets simultaneously on 
different tracks displaying all targets on one heat map, which would allow JIATF South planners 
to analyze more complex target scenarios. Figure 3.4 shows a scenario with four targets on the 
water en route to various destinations in Mexico and Central America. Figure 3.4a depicts the 
four targets 15 hours into the scenario and Figure 3.4b shows three targets, as one of the four has 
reached it’s destination. The introduction of departure time uncertainty is evident in the long 
blue channel in Figure 3.4b. This particular target has a departure window of 48 hours thereby 
spreading the probability mass across the entire voyage. We refer to these cases as "channel" 
cases. The optimization model handles these cases differently and requires additional inputs. 
We currently exclude these cases from analysis. Following the construction of a more general 
probability model, we begin the process of analyzing the probability maps to reproduce the 
target voyage characteristics. 
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Figure 3.4: A multiple target scenario produced in the enhanced probability model. The more 
general probability model is able to consider several targets simultaneously, allowing JIATF 
South planners to analyze more complex target scenarios. The long "channel" is a special case 
when the departure window is large. The optimization model handles these differently and we 
exclude these types of cases from analysis. 
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CHAPTER 4: 
Path-finding Model 


A strong focus of this research analyzes how to create the target’s track from the probability map 
generated by the enhanced probability model. Specifically, we determine the target’s earliest 
and latest departure times, earliest and latest arrival times, speed on each leg, course on each leg, 
location of waypoints, time of change in course or speed. Area of Uncertainty (AOU), departure 
location, and arrival location. This data is then fed into the optimization model described in 
Section 2.3 to determine optimum employment of JIATF South assets to maximize the expected 
amount of drugs seized. Additionally, the analytic probability model provides the ability to 
verify more complex simulation models such as An et al. (2012)’s extension of piracy model to 
counter drug operations. The path-finding model assumes that the probability data is the only 
data available where p{x,y,t) is the probability that a target is located in a grid cell centered at 
(jc,y) at time t. From this data, we derive all other parameters for the optimization model. 

4.1 Optimization Model 

The Pietz and Royset (2013) optimization model dictates the output required from our path¬ 
finding model. This optimization model is a path-constrained optimal search problem in con¬ 
tinuous space and time (Pietz & Royset, 2013). The model determines the optimal routing of 
search assets in the area of operations to maximize the expected drug load. The area of oper¬ 
ations is modeled as a transportation network, in which the nodes are search regions and the 
arcs are the searcher’s path between moving search regions. These search regions are areas of 
uncertainty concerning the target’s location, thus knowledge of target parameters is necessary 
to run the optimization model. The inputs associated with the targets necessary to run the opti¬ 
mization model are listed in Table 4.1 and they correspond to those defined in Pietz and Royset 
(2013). 

The optimization model assumes that the targets travel in a piecewise linear motion at a known 
constant speed (Pietz & Royset, 2013). The path-finding model derives the uncertainty data 
listed in Table 4.1 for input into the optimization model. Initially we identify any waypoints, 
which correspond to changes in a target’s course or speed and divide the path into the appro¬ 
priate number of legs. Subsequent analysis determines velocity, initial position, final position, 
expected departure time, departure time uncertainty, departure/arrival position uncertainty, and 
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Expected departure time of target j Xj 

Time uncertainty range of target j Xj 

Expected departure location of target j pj 

Expected arrival location of target j pj 

Departure/arrival location uncertainty range of target j pj 
Expected value of detecting target j qj 

Speed of target j Uj 

Table 4.1: The inputs required to run the optimization model creating a search vehicle allocation 

plan. (Erom Pietz and Royset, 2013). 


AOU for each leg of the target’s path. It is possible some of the uncertainty data in Table 4.1 
is directly available from JIATE South databases or intelligence. However, in a worst case sce¬ 
nario, planners have only the probability map and must determine the aforementioned inputs 
solely from this map to produce a search vehicle allocation plan via the optimization model. 

4.2 General Path-finding Model Approach 

The path-finding model approximates the path taken by the target as a sequence of legs where 
the collection of all legs produces the target’s entire track from departure location to arrival 
location. We define a leg as a segment of the target’s track that has a constant heading, speed, 
and AOU. This creates a piecewise linear path for the target. Points in time in which the target’s 
speed, heading, or AOU changes signal the ending of one leg and the beginning of the next. We 
address new legs based on changes in the AOU once all the constant speed-heading legs are 
identified. Eigure 4.1 shows a potential path from departure location to arrival location for one 
target. In this case, the target’s path consists of two legs in which the target changes heading 
at the waypoint. The mechanics of the path-finding model automatically identifies where and 
when a target changes heading or speed. This allows us to generate the piecewise linear path 
required in the optimization model. 

On each leg we must produce linear tracks with a constant speed, heading, and AOU. We deter¬ 
mine these inputs for each target separately and compile the inputs by leg for the optimization 
model. We apply regression techniques to generate piecewise linear paths through the proba¬ 
bility mass. Specifically, the path-finding model produces two linear regression equations of 
the target’s location, one for longitude and one for latitude. The coefficients associated with the 
two linear regression equations are then used to calculate the required quantities. The associ¬ 
ated mathematical equations are discussed further in Section 4.4 and numerical illustrations are 
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Figure 4.1: Illustration of one target traveling along two legs from departure to arrival location. 
This track consists of two legs with a heading change at the waypoint. 


provided in Sections 4.5 and 4.6. 

From the probability model, we analyze a set of non-zero probabilities that are spread over a 
cellular grid in the AOR. For each time period, the probabilities correspond to the center of a 
grid cell by their associated latitude and longitude. We define n as the number of grid cells in 
the AOR, pi{t) as the probability the target is in cell i at time t, and (A:,,y,) as the center point 
of grid cell i. Figure 3.2 illustrates four different hours of the probability mass in the AOR with 
each grid cell colored based on the target’s probability of being in that cell at the associated 
time. This probability mass changes size and shape and moves over time as the target travels 
from the departure location to the arrival location. At each time step, we calculate the expected 
value of the target’s longitude or x component where £ [A] = Pi{t)xi and the expected value 
of the target’s latitude or y component where £’[7] = YJi=i Pi{t)yi- 

Since the probability mass is spread over several grid cells during each time step, weighted 
target tracks could be produced but for simplicity and limitations of the techniques employed, 
we use expected values. 

The resultant expected x and y coordinates are then plotted in Figure 4.2 for a visual picture of 
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the target’s expected track. Figure 4.2 illustrates the same track from Figure 4.1 but uses the 
expected values of X and Y. As in Figure 4.1, Figure 4.2 shows two legs with a heading change 
near 86°W 3°N (274°, 3°). 


Expected Target Track 



Figure 4.2: Plot of one target’s expected track with two legs in the AOR. 

Based on the plot of latitude versus longitude, we can easily see a change in heading but dis¬ 
cerning departure and arrival locations or speed changes is not possible. The necessary time 
element is not depicted in Figure 4.2. This leads us to plot longitude versus time and latitude 
versus time as shown in Figures 4.3 and 4.4 respectively. These plots become the primary means 
of determining the necessary target parameters in Table 4.1. 

Our regression analysis utilizes the probability data between the latest departure time and the 
earliest arrival time, which corresponds to the time that the entire probability mass is on the 
water. These points are specifically highlighted in Figures 4.3 and 4.4. Any change in slope 
between these points signifies a change in the target’s speed, heading, or possibly both. Changes 
in slope outside this interval may produce spurious changes in course or speed due to endpoint 
effects and are not important for optimization. Figures 4.3 and 4.4 both show a change in slope 
around hour 45, signifying a change in the target’s motion. The points where a change in slope 
exists in either the x direction or y direction indicates where one leg ends and the next leg 
begins. 

Detecting two legs is obvious in Figure 4.1 because there is a considerable change in course. 
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Expected Target Track in Degrees Longitude Over Time 



Figure 4.3: Plot of one target’s expeeted traek with two legs in degrees longitude over time. 


Expected Target Track in Degrees Latitude Over Time 



Figure 4.4: Plot of one target’s expeeted traek with two legs in degrees latitude over time. 


However, a ehange in speed with a eonstant eourse ean not be determined with a statie figure. 
Consider the example in Figure 4.5 where the target is traveling in a straight line from the 
departure loeation to the arrival loeation. 

Based on Figure 4.5, it appears one leg exists in this target’s path. Looking at the expeeted traek 


29 









Figure 4.5: Illustration of a potential target traveling in a straight line path from the departure 
location to the arrival location. 

in Figure 4.6 we are inclined to the same conclusion. 

Expected Target Track 
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Figure 4.6: Plot of the one target straight line path in degrees latitude versus degrees longitude. 

A closer look at this figure shows the spacing between expected locations increasing. This 
increase in space between expected locations signifies the target has changed speed. We are still 
missing the time element and cannot tell whether the target is traveling northwest or southeast 
but we can tell speed has either increased or decreased during the journey. When we consider 
the plots in Figures 4.7 and 4.8, it is evident the target is traveling northwest and increases speed 
during the journey based on the increase in slope at hour 25. 
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Expected Target Track in Degrees Longitude Over Time 
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Figure 4.7: Plot of the one target straight line path in degrees longitude over time. 


Expected Target Track in Degrees Latitude Over Time 
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Figure 4.8: Plot of the one target straight line path in degrees latitude over time. 


This ease elearly illustrates the neeessity of analyzing the x and y eomponents of the target’s 
motion separately and the need to employ a teehnique eapable of produeing pieeewise linear 
equations. 
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4.3 Path-finding Model 


We desire to provide JIATF South with a path-finding model eapable of performing this analysis 
automatieally. To meet this need we employ a teehnique known as MARS. MARS, as deseribed 
in Seetion 2.4 automatieally fits a pieeewise linear regression model to the probability data by 
automatieally seleeting "knots" or points of ehange in eourse or speed. 

MARS uses a forward model building and baekward term deletion proeedure to produee the 
final model (Hastie et ah, 2009). During the forward model building proeedure, MARS uses a 
set of basis funetions as shown in Equations 4.1 and 4.2 as input. 


{x-t) + 


x — t, 

0 , 


ifx>t, 

otherwise 


(4.1) 


{t-x)+ = 



ifx<t, 

otherwise. 


(4.2) 


For eaeh input, MARS forms a set of refleeted pairs, {Xj —1)+ and (t — Xj)^, for eaeh value of 
that input (Hastie et ah, 2009). The resulting model is of the form 


M N 

/(A)=j8o+ £i8„(A-t^)++£j8„(t„-A)+ (4.3) 

m=l M=1 

where tm and are the eandidate times of eourse or speed ehanges. Additionally, the model 
form ean inelude a product of two or more basis functions (Hastie et ah, 2009). For a particular 
basis function, MARS estimates the associated coefficient, /3, by minimizing the residual sum- 
of-squares (Hastie et ah, 2009). The core functionality of MARS lies in the selection of the 
basis functions. Each step through the forward model building process a new basis function 
pair and its products of all functions in the model are evaluated for inclusion in the model. The 
term that produces the largest decrease in error is then added to the model. For a full description 
of MARS see Friedman (1991) and Hastie et al. (2009). 

Consider the target scenarios presented in Figures 4.1 and 4.5. As previously discussed, these 
targets follow tracks that consists of two legs identified by a course change in the expected track 
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and a speed change respectively. Utilizing MARS, the path-finding model automatically locates 
changes in slope as seen in Figures 4.9 through 4.11 and produces final fitted models of the form 


^ — Po-\- — 45)+ -|- ^ 2(45 —1)+ 


(4.4) 


y = j8o + /3i(t-37)++/32(t-49)+ + j 83 ( 49 -t)+. (4.5) 


Expected Track in Degrees Longitude 
Fith with MARS and Linear Regression 



Figure 4.9: Plot of one target’s track with two legs in degrees longitude over time. The resulting 
MARS model fits the path extremely well, whereas ordinary least squares regression is unable 
to automatically identify the course change. 

Figures 4.9 through 4.11 show the resulting MARS model fit and illustrate the poor resulting 
fit by using standard linear regression. Standard regression does not adequately identify the 
slope changes, thus it does not produce the piecewise linear track that MARS does. Note, 
Figure 4.10 displays two slope changes and Equation 4.5 includes different basis functions than 
Equation 4.4. In fitting models in the x and y directions separately, MARS produces a best fit 
model irregardless of the number of basis functions in the other direction. A few reasons exist 
for the additional knot and are discussed further in Section 4.6 as well as the case parameters 
and associated model output. 
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Expected Track in Degrees Latitude 
Fit with MARS and Linear Regression 



Figure 4.10: Plot of one target’s track with two legs in degrees latitude over time. The resulting 
MARS model fits the path extremely well, whereas ordinary least squares regression is unable 
to automatically identify the course change. 


Expected Target Path 
Fit with MARS and Linear Regression 



Figure 4.11: Plot of one target’s track with two legs in degrees latitude versus degrees longi¬ 
tude. The resulting MARS model fits the path extremely well, whereas ordinary least squares 
regression is unable to automatically identify the course change. 
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At this point, the basis functions represent potential waypoints along the target’s route. Next, 
the path-finding model checks the validity of these waypoints by first evaluating the change 
in slope of reflected pairs, and then evaluating the time between waypoints as compared to 
the time required for the probability mass to complete the associated turn. These concepts are 
demonstrated in greater detail in Sections 4.5 and 4.6 

Once we know which waypoints are valid, we construct a new data set to account for the way- 
points. We add a new data field, wpi, for each valid waypoint where wpi equals the value 
of the associated basis function. If t^p is the time of the respective waypoint, then we define 

wp = 

The path-finding model uses this data set and ordinary least squares regression to produce the 
final model, which is of the form 


X = 


j 

Po+pit + Y^PiWPi 


i={ 


(4.6) 


and 


y = Po + pit + Y^Piwpi (4.7) 

i=l 

where j represents the number of waypoints. The additional data fields, identified by wpi, enable 
ordinary least squares regression to produce a piecewise linear fit. 


4.4 Input Generation 

The probability data and the path-finding model produce the uncertainty parameters listed in 
Table 4.1 to generate a search vehicle allocation plan. Recall we have uncertainty in departure 
location, arrival location, departure time, and velocity and in a worst case scenario JIATF South 
planners will only have the probability map to generate the needed parameters for optimal 
allocation of their search assets. The following section describes how each of the inputs are 
generated from either the probabilities or the model output. The desired end-state equates to a 
list of parameters for each target and leg of the associated target’s path. 
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4.4.1 Times 

The introduction of T, the target’s random departure time, in the enhanced probability model 
produces uncertainties in the departure time and subsequently the arrival time for a target j. The 
uncertainty results in a range of possible departure times between [7/,r„] and a corresponding 
range of arrival times. The specific times for each leg required in the Pietz and Royset (2013) 
optimization model are listed in Table 4.2. Additionally, we define the earliest departure time 
of target j from the origin as and the latest arrival time at the destination as 

for computational purposes. 


Expected departure time of target j Xj 
Time uncertainty range of target j Xj 
Latest departure time of target j 
Earliest arrival time of target j 

Table 4.2: Times required on each leg for the optimization model to generate a search vehicle 
allocation plan. (Erom Pietz and Royset, 2013) 


Pietz and Royset (2013) discuss calculating Tf'" and from the uncertainty data, T,- and f 
Our path-finding model first determines ^ ^max^ ^latestarr 

last leg respectively and then we calculate Xj and Xj. It is important to note that these initial 
time determinations are for the first leg and the last leg and are determined directly from the 
probability data. Each leg of the target’s path requires the identification of the times listed in 
Table 4.2. MARS generates these times for the intermediate waypoints, whereas the probability 
data generates the times associated with the origin and destination. Eigure 4.12 illustrates the 
associated times. In this case more than one leg exists, thus the intermediate departure and 
arrival times are determined through regression analysis. 

The earliest departure time is the first time period that any non-zero probability exists whereas 
the latest departure time is the first time period that the sum of the probabilities is equal to 1. 
These are expressed in Equations 4.8 and 4.9 respectively. 


earliestdep 

Xj = mrn • 


i=l 


^min 


= min 


= i| 


(4.8) 

(4.9) 
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Expected Target Track in Degrees Longitude Over Time 



Figure 4.12: The initial time determinations in the path-finding model represent the earliest and 
latest departure times on the first leg and the earliest and latest arrival times on the last leg. 


Similarly for the arrival times, the earliest arrival time is the maximum time that the probabilities 
sum to 1 and the latest arrival time is the maximum time the probabilities sum to greater than 
zero. These are expressed in Equations 4.10 and 4.11 respectively. 


T;“ = max|t|£p,(t) 



(4.10) 


-latestarr 


= max 


tlY.Piit) 

i=l 



(4.11) 


It is helpful to view the origin as a source and the destination as a sink. The probability mass 
begins to grow from the source and eventually 100% of the probability mass has left the source, 
thus representing the earliest and latest departure times respectively. As the target travels across 
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the AOR, the probability mass begins to disappear into the sink until all of the mass has finally 
collected in the sink, thus representing the earliest and latest arrival times. When only one leg 
exists for a target, all times are determined directly from the probability data. The departure 
time uncertainty range of target j, Xj, is simply the difference between the earliest departure 
time and latest departure time. Based on ^nd t”'”, we determine fj and Xj using 

Equations 4.12 and 4.13 respectively. 


~ _ -min _ earliestdep 


(4.12) 




mm _ earliestdep\ 

j V ^ 


(4.13) 


If the departure window is large enough, the probability mass can grow from the source to the 
sink without ever fully leaving the source. These cases are referred to as "channel" cases and 
the optimization model handles them slightly different. An example is shown in Figure 3.4b. 
At this point, we only focus on cases where the probability mass develops entirely on the water 
before beginning to absorb into the sink. For purposes of the path-finding model, we are mainly 
concerned with the latest departure time and the earliest arrival time. In other words, the path¬ 
finding model only uses the probability data when the probabilities sum to 1, which is the 
probability data between the first T™'” and the last t™". 

4.4.2 Position 

Each position is composed of a horizontal component and a vertical component and the random 
variable D introduces uncertainty in the departure location on the interval The regres¬ 

sion models, when run between the first t™'” and the last result in a set of fitted values 
between these times. Essentially, they are an estimate of the target’s horizontal and vertical 
position in the AOR for each time step between departure and arrival. The optimization model 
requires position inputs for each leg, specifically a starting position and an ending position. 
Pietz and Royset (2013) define the origin of the AOR as 0° N (0°) and 100° W (260°) where 1° 
in any direction equals 60 nm. All positions are measured from this origin. 


The fitted values from the regression model directly correspond to the position of target j. 
Knowing T™'" and for each leg allows us to determine the target’s position by substi¬ 
tuting and x1^^ into the final models. The resulting {x ^min , 3^ ^min ) and are 
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then converted to nautical miles from the AOR origin where pj{x) = (xi—xq)* 60 and pj(y) = 
iyt-yo)*60. 

Pietz and Royset (2013) define pj and pj as the expected departure location and the expected 
arrival location of target j. In other words, pj and pj are the target’s position at Tj and the last 
^max j-espectively. An accurate estimate of the expected departure location is required for the 
optimization model, but an accurate estimate of the expected arrival location is not required 
as long as its along the same track. Since we regress on the probability data between the first 
^mm ^max^ estimate, based solely on regression, of p, is always further from the 

true origin. This fact does not impact the path-finding model, however; to correct this, we use 
the vertical and horizontal speeds to move the position at the first T™'" the distance of one-half 
the departure time uncertainty times the respective speed. For example, with a departure time 
uncertainty of 4 hours and a vertical speed of 10 knots and a horizontal speed of 12 knots, the 
pj for the first leg moves 20 nm in the vertical direction and 24 nm in the horizontal direction 
towards the true origin. 


4.4.3 Speed 

The optimization model assumes the searcher knows the speed of target j on each leg of the 
path and defines this speed as Uj (Pietz & Royset, 2013). This assumption then applies to 
the entire leg without considering the possibility of speed changes during the leg. Recall, a 
change in speed or heading is cause for the introduction of a new leg since the optimization 
assumes a constant speed on each leg. Considering the worst case scenario, where only the 
probability data exists, we desire to determine the target’s speed on each leg. The coefficients 
of the regression model in Equations 4.6 and 4.7 are used to calculate speed. The coefficients 
in the linear models, excluding the intercept (jSq), represent the speed in the x and y directions. 
Overall target speed is then defined as Uj = 


In a multiple leg case, the regression techniques produce a piecewise linear model for the entire 
journey where, depending on the time, certain coefficients are active. Recall the MARS dis¬ 
cussion in Section 4.3, the basis functions equal zero when they evaluate to a negative quantity. 
Therefore, for times when a basis function equals zero the respective coefficient does not con¬ 
tribute to the target’s speed. The coefficients only contribute to the speed calculation when the 
basis function evaluates to a non-negative number. Knowing the T™'" and for each leg, the 
path-finding model determines which coefficients are active and based on the active coefficients 
computes the speed for the respective leg. 
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4.4.4 Area of Uncertainty 


An AOU results from uncertainty in both the departure time and the departure location. This 
area is essentially the probability mass moving across the water. Just as a constant speed on each 
leg is required, so is a constant area along each leg. Changes in the AOU also would indicate 
the need for additional legs. Determining the area of the probability mass as it moves across 
the water is rather straightforward. We can simply count the number of grid cells that contain 
probabilities greater than 0. Knowing the size of each grid cell and the number of grid cells 
with a non-zero probability, we simply multiply the two and determine the associated AOU. 
This method while straight forward does not result in a constant area over time. In general, the 
area changes every time step and the user may not require 100% coverage of the AOU. We 
explore several methods to determine the AOU. 

The first approach determines the area of the minimum bounding box at each time step of 
the target’s path. This method captures the entire probability mass and results in the largest 
AOU. Additionally, it is possible for the minimum bounding box to capture locations with zero 
probability, as evident in Figure 4.13. Therefore, the second approach determines the convex 
hull of the probability mass. Again, the convex hull captures all of the probability mass but with 
minimal area. Figure 4.13 depicts the minimum bounding box and convex hull at several time 
steps. These methods successfully bound the entire probability mass and capture 100% of the 
area. 

The third approach uses a greedy algorithm to capture a user defined amount of the probability. 
Throughout this thesis we assume a probability threshold of 90%. For each time between the 
first and the last the greedy algorithm selects the node or cell with the highest prob¬ 
ability and then creates a set of cells that neighbor the initial cell. The algorithm proceeds by 
examining the set of neighbors and selects the neighbor with the highest probability and adds 
it to the AOU set. Next it calculates the cumulative probability of the AOU set and proceeds to 
look at the set of neighbors including any new neighbors, again adding the neighbor with the 
highest probability. This procedure continues until the cumulative probability is greater than or 
equal to the selected threshold. Once we reach the specified threshold, we simply determine the 
amount of cells in the set and multiply by the area of a single cell to calculate the total AOU. 
Figure 4.13 shows the grid cells that are included in the greedy algorithm’s AOU set. 

Figure 4.14 shows the results of the three methods over time. The AOU found using the mini¬ 
mum bounding box tends to result in a considerably higher area since it encompasses the entire 
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(a) Hour 20 



(b) Hour 28 
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Figure 4.13: Examples of the minimum bounding box, eonvex hull, and greedy algorithm 
applied to a target’s probability mass over several time steps of a target’s path. The minimum 
bounding box and eonvex hull eneompass the entire probability mass at each time step, while 
the greedy algorithm results in the smallest AOU. 
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probability mass and additional area with zero probability. The convex hull performs slightly 
better but it is still greater than the greedy approach. By taking a probability threshold of 90%, 
we are able to significantly reduce the size of the AOU while still ensuring we capture a signif¬ 
icant amount of the probability mass. Additionally, the greedy algorithm allows for flexibility 
and if desired, the user may select a threshold of 100%. 


Target Area of Uncertainty 



Figure 4.14: Comparison of the AOUs resulting from the three methods explored. The greedy 
algorithm consistently produces areas less than the other two approaches, while the minimum 
bounding box results in significantly higher areas. 


With the AOU specified for each time step of the target’s path, we now desire to create a piece- 
wise constant area for the entire journey. The area analysis utilizes the functionality of MARS 
and automatically fits a piecewise linear model to the area data. The model smooths the area 
data but does not produce the piecewise constant areas required. Figure 4.15 shows the MARS 
model fit to the areas from the greedy algorithm. This implementation of MARS forces a mini¬ 
mum span between knots equal to the departure time uncertainty. The minimum span prevents 
MARS from introducing knots at points just before or after a new leg. 

Additionally, for each constant speed and heading leg, we build a separate model since we 
desire to create the minimum number of new legs based on area, and we know the number of 
legs based solely on speed and heading changes. Figure 4.15 depicts the MARS models fit to 
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MARS Model Fit to Greedy AOU 



Time (Hours) 

Figure 4.15: A MARS model smooths the AOU over time. A different model is built for each 
leg with a constant speed and heading. The fitted values from the model are further analyzed to 
produce a piecewise constant AOU. 


each leg of the target’s path. To determine the number of additional legs based on a change in the 
AOU, we calculate the maximum number of legs that we can fit in each constant speed-heading 
leg based on two different criteria. First, we use the minimum span to determine the number of 
additional legs that can be fit into each constant speed-heading leg. This is simply the length 
of the leg in time divide by the departure time uncertainty. The second method takes a user 
selected percent change in area, which is currently set at 25%, and determines the number of 
times a 25% change in area occurs. Then we take the minimum between the two methods. This 
gives us the number of new legs to insert based on area changes. We perform this procedure 
on each constant speed-heading leg. Considering the slope of the MARS model, we assign 
areas to the new legs at the maximum area of the new interval. Figure 4.16 shows the resulting 
piecewise constant area. 

In this particular case, there are two constant speed-heading legs where the course changes at 
hour 30. Based on the AOU analysis, we divide the first constant speed-heading leg into two 
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Piecewise Constant AOU 
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Figure 4.16: Illustration of a piecewise constant AOU. 


segments. The first segment ends at hour 21, at which point, the AOU increases. The second 
constant speed-heading leg remains unchanged because the AOU does not change by more than 
25%. 


4.5 Straight Line Model and Illustration 

Assume there is one target of interest that travels from the origin to the destination with a 
constant speed and heading. Figure 4.17 shows the target’s track from departure location to 
arrival location. We illustrate the path-finding model by first entering the case parameters into 
the enhanced probability model. The case parameters as shown in Table 4.3 include uncertainty 
in departure and arrival locations and departure time. 

Running the probability model produces the probability data over time necessary for the path¬ 
finding model. Figure 4.18 depicts several hours of the target’s location as it moves through the 
AOR. In a worst case, this is all the JIATF South planners have and they desire to determine the 
information in Table 4.3. 
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Figure 4.17: Illustration of a target’s traek from departure to arrival loeation. The target follows 
a direet path at a eonstant speed to illustrate the simplest path-finding ease. 


Approximate payload 

Speed 

Origin 

Departure uneertainty 
Destination 
Arrival uneertainty 
Departure time 
Departure time uneertainty 
Intel eertainty 


5,000 kg 
15 knots 
76.3° W 9° N 
±30 nm, Uniform 
83.5° W 14° N 
±30 nm, Uniform 
33 hours 

±4 hours. Uniform 
95% 


Table 4.3: Parameters for a simple straight path from origin to destination with eonstant veloeity. 
Uneertainty data is listed for departure and arrival loeations and departure time. 


In this ease, the probability model generates the heat map for 78 hours using the assoeiated 
parameters from Table 4.3. In the eontext of this thesis, hour 1 represents the first simulated 
time step and hour 78 represents the last time step. This results in 78 probability data files, one 
for eaeh time step, that we read into the path-finding model. Eaeh file eontains a probability 
for eaeh grid eell in the AOR eorresponding to the likelihood the target is in that grid eell at the 
assoeiated time. The path-finding model analyzes the probabilities to establish the times listed 
in Table 4.2 using Equations 4.8 through 4.13. Table 4.4 lists the resulting times. 
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(a) Hour 37 


(b) Hour 52 



(c) Hour 64 


Figure 4.18: A simple one target seenario produeed in the enhaneed probability model. The 
snapshots are at the latest departure time (hour 37), an intermediate point (hour 52), and at the 
earliest arrival time (hour 64). 



Once the first departure time (t™'^) and last arrival time (t™“) are known, the path-finding 
model proceeds by applying MARS to the probability data between these times. The resulting 
form of the MARS model is 


= 280.42-0.205(1-49)++ 0.205(49-0+ (4.14) 


and 
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Time 

Actual 

Model Result 

^earliestdep 

29 hours 

29 hours 

First 

37 hours 

37 hours 

Fast 

64 hours 

64 hours 

^latestarr 

NA 

72 hours 

f 

8 hours 

8 hours 

T 

33 hours 

33 hours 


Table 4.4: Case time information from path-finding model as compared to the actual inputs for 
a target traveling along a straight path. 


yt = 12.13+ 0.142(t-55)+-0.142(55-t) + . (4.15) 

From Equations 4.14 and 4.15, we see MARS introduces the basis functions (t —49)+ and (49 — 
t)+ in the x direction and (t — 55)+ and (55 —1)+ in the y direction. Further examination reveals 
the coefficients associated with the reflected pairs have equal magnitude but opposite sign. This 
is a result of the MARS model fitting procedure. Continuing the path-finding algorithm with 
this model results in a track composed of three legs but the same speed and heading on each leg. 
In this particular case, it is easy to determine that these points are not truly changes in speed or 
heading. In a more general case, it is necessary to determine if a point introduced by MARS is 
really a change in speed or heading. 

We apply a filter to the knots or cuts introduced by MARS. The filter first examines the change 
in slope at each reflected pair. If the absolute difference in slope for both models changes by 
less than 0.01 the filter does not add them to the expected track and path-finding continues with 
the expected target locations as defined in Section 4.2. In cases where the reflected pairs or 
cuts are legitimate, we add an additional data field to the expected track data. This data field 
is identified by the associated time and populated with the result of the basis function for each 
time step in the path. For the moment let us consider the reflected pair at hour 49 as a legitimate 
time to split the journey in the x direction into two legs. The expected track data would then 
get an additional data field containing the result of (t — 49)+ and the final regression model is 
constructed with the new data. 

In the current case, the filter determines that the points at time 49 and 55 are not valid based on 
the change in slope. These are not added to the expected data and the algorithm continues the 
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path-finding procedure using the expeeted loeations in the x and y direetions. At this point of 
the proeedure, we know the target is traveling from the departure loeation to the arrival loeation 
on a straight path at a eonstant speed. Reeall that the optimization model requires a eonstant 
AOU so we proeeed to determine the AOU throughout the target’s path. 

The associated plot of the AOU and the resultant MARS model fit is shown in Fiugre 4.19. The 
overall ehange in area during the leg is less than the seleeted threshold of 25%, thus the AOU 
is chosen as the maximum area during the path. In this case, we capture 90% or more of the 
probability mass during the entire journey. 


Piecewise Constant AOU 



40 45 50 55 60 65 


Time (Hours) 


Figure 4.19: Resulting AOU throughout one target’s straight path and assoeiated MARS model 
fit to the greedy area data. 


With the entire traek examined for speed, heading, and area ehanges, we proeeed to build the 
final ordinary least squares model using the expected location and time data. For this particular 
case with no waypoints and no ehanges in the AOU, the final model takes the following form 


T = 290.5-0.205t 


(4.16) 
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and 


3; = 4.3+ 0.142?. 


(4.17) 


Based on the times in Table 4.4 and Equations 4.16 and 4.17, we determine the remaining 
parameters utilizing the equations described in Section 4.4. Table 4.5 compares the results of 
the path-finding model with the associated inputs to the probability model. The path-finding 
model determines the target’s departure location, speed, and associated departure and arrival 
times extremely well. When comparing the destinations, remember that the path-finding model 
uses the probability data for the time periods when the probabilities sum to one. Thus, the model 
determines an arrival location that is actually on the water vice on the coast as depicted by the 
actual location. An accurate estimate of the origin is required for the optimization model, but 
an accurate estimate of the destination is not required as long as its along the same track. 


Parameter 

Actual 

Model Result 

Speed 

15 knots 

14.981 knots 

Origin 

76.3° W 9° N 

76.3° W 9° N 

Destination 

83.5° W 14° N 

82.657° W 13.414° N 

AOU 

7200 nrn^ 

7582.65 nni^ 

^earliestdep 

29 hours 

29 hours 

^min 

37 hours 

37 hours 

^max 

64 hours 

64 hours 

f 

±4 hours 

±4 hours 

T 

33 hours 

33 hours 


Table 4.5: Comparison of actual parameters and path-finding model results for the simple target 
scenario. 


Figure 4.20 shows the final model fit to plots of the expected track through the AOR and over 
time. 


4.6 Multiple Leg Model and Illustration 

Recall the target with two legs from Figure 4.1. Table 4.6 lists the associated parameters that 
we enter into the enhanced path-finding model. 

Again we show in Figure 4.21 several hours of the probability mass as it moves through the 
AOR. The probability model generates a heat map for 120 hours of the target’s journey and the 
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Final Model Fit to Expected Location for 
One Target on a Single Leg 




Case 1 



Figure 4.20: Output of the path-finding model, showing the target’s track through the AOR and 
the target’s track in both the x and y direction with respect to time. 


resulting probability files feed into the path-finding model. Based on the probability data, we 
determine the initial and final departure and arrival times which are listed in Table 4.7. 
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Approximate payload 

Speed 

Origin 

Departure uncertainty 
Waypoint 

Waypoint uncertainty 
Destination 
Arrival uncertainty 
Departure time 
Departure time uncertainty 
Intel certainty 


2,500 kg 
15 knots 
IT W 5° N 
±30 nm, Uniform 
86° W 3° N 
±60 nm 

92.2° W 14.5° N 
±50 nm, Uniform 
10 hours 

±3 hours. Uniform 
50% 


Table 4.6: Parameters for a multiple leg path from origin to destination. Uncertainty data and 
associated distributions are listed for departure location, arrival location, and departure time. 


Time 

Actual 

Model Result 

^earliestdep 

7 hours 

7 hours 

First t"“" 

13 hours 

13 hours 

Last 

96 hours 

91 hours 

^latestarr 

NA 

107 hours 

f 

6 hours 

6 hours 

T 

10 hours 

10 hours 


Table 4.7: Case time information from path-finding model as compared to the actual inputs for 
a target traveling along a path with multiple legs. 


Next, we proceed with the application of MARS between time 13 and time 91. The resulting 
form of the MARS model is 


f, = 274.36-0.123(t-45)+±0.249(45-0+ (4.18) 


and 


ft = 2.77 ± 0.052(t - 37)+ ± 0.170(t - 49)+ ± 0.058(49 - 0+- (4.19) 

In this case, MARS introduces waypoints at hours 37, 45, and 49. We examine the validity of 
each waypoint by looking at the change in slope and the maximum time for the entire probability 
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- 0.0 

I 


Figure 4.21: A one target, multiple leg scenario produced in the enhanced probability model. 
Snapshots are from various times during the target’s route. 


mass to travel through the respective waypoint. This time is then compared to the amount of 
time until the next waypoint. If the maximum time is greater than the time to the next waypoint, 
then we discard the next waypoint and assume it is a artifact of the probability mass moving 
through a turn. 

Here, the change in slope is significant, thus solidifying that we have at least one valid waypoint. 
Next we determine which of the three are valid. We assume that the first waypoint, chronologi¬ 
cally, is valid and check if the maximum time to move through the turn is less than the time to 
the next waypoint. Figure 4.22 illustrates the probability mass at hour 45. The current waypoint 
(at hour 45) and the next waypoint at hour 49 are highlighted. We desire to determine the length 
of time for this probability mass to move through the hour 45 waypoint. 
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Probability Mass at Hour 45 



Figure 4.22: Depiction of the probability mass in the one target multiple leg case. The prob¬ 
ability mass is shown at hour 45, which is the time of the current waypoint. We rotate this 
probability mass to align with the target’s course en route to the next potential waypoint. 


Specifically, we rotate the probability mass along the target’s path to determine the maximum 
distance between x coordinates. This rotation aligns the probability mass with the target’s track, 
allowing us to use the target’s overall speed in the time calculation. This distance is then divided 
by target speed to produce the time to complete the turn. At hour 37, the time to complete the 
turn is 6.8 hours, which is less than the time to the next waypoint (8 hours), thus the waypoint 
at hour 45 is kept. The turn time at hour 45 is 7.4 hours and is greater than the time to the next 
waypoint (4 hours), therefore we discard hour 49 as a waypoint and conclude the target path 
consists of three legs based on speed and heading changes. Figure 4.23 illustrates the rotated 
probability mass and the associated horizontal distance. The final path based on constant speed 
and heading legs contains two valid waypoints with the first at hour 37 and the second at hour 
45. 

The AOU analysis follows which determines if the path requires additional legs based on sig¬ 
nificant changes in the AOU. The associated plot of the AOU and the resultant MARS model 
fit is shown in Figure 4.24. In this case, we analyze each leg separately using MARS. 

The resultant MARS models fit three different linear piecewise segments to the area data. The 
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Rotated Probability Mass at Hour 45 
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Figure 4.23: Depiction of the probability mass in the one target multiple leg case. The rotated 
probability mass is shown at hour 45, which is the time of the current waypoint. The rotated 
probability mass aligns with the target’s course, which allows us to calculate the time for the 
target to move the maximum x distance using overall target speed. 


first segment begins at hour 13 and ends at hour 37. During this time, the area changes by more 
than 25%, thus the model divides the constant speed heading leg into two legs based on the 
change in area. The area on the other two legs, from 37 to 45 and from 45 to 91, changes by 
less than 25% and are not further divided. Figure 4.25 depicts the associated times and areas of 
the new legs. 

With the addition of new legs based on area changes, the path now consists of four legs. The 
new leg does not change the speed or heading of the target but only represents the location that 
the AOU increases by more than 25%. As a result, the final model accounts for the speed and 
heading changes identified by MARS but it is unnecessary to insert the points of area change. 
The final path-finding model is of the form 

xt = 285.4 - 0.242t - 0.024(t - 37)+ + 0.145(t - 45)+ (4.20) 


and 
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MARS Model Fit to Greedy AOU 



Figure 4.24: Resulting AOU throughout one target’s multiple leg path and associated model fit 
to the greedy area data. 


Piecewise Constant AOU 



Figure 4.25: Resulting piecewise constant area for one target multiple leg case. 
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= 5.47 - 0 . 050 ? - 0 . 055 (? - 37 )+ + 0 . 320 (? - 45 )+. 


( 4 . 21 ) 


The mechanics to determine the remaining case parameters are the same as those in Section 4.5 
and Table 4.8 lists the model output and the actual parameters. Figure 4.26 shows the final 
regression model fit to the expected track over time. 


Parameters 


Leg 1 


Leg 2 


Actual 

Model 

Actual 

Model 


Speed 

15 knots 

14.843 knots 

NA 

14.843 knots 


Origin 

IT W 5° N 

77.015° W 4.964° N 

NA 

80.649° W 4.211° 

N 

Destination 

86° W 3° N 

80.649° W 4.211° N 

NA 

83.555° W 3.608° 

N 

AOU 

5400 

5697 

NA 

7422 


^min 

13 

13 

NA 

25 


^max 

46.88 

25 

NA 

37 


f 

6 

6 

NA 

6 


T 

10 

10 

NA 

25 




Leg 3 


Leg 4 


Parameters 

Actual 

Model 

Actual 

Model 


Speed 

NA 

17.175 knots 

15 knots 

14.754 knots 


Origin 

NA 

83.555° W 3.608° N 

86° W 3° N 

85.684° W 2.763° 

N 

Destination 

NA 

85.684° W 2.763° N 

92.2° W 14.5° N 

91.264° W 12.603'^ 

’N 

AOU 

NA 

9184 

10800 

10077 


^min 

NA 

37 

46.878 

45 


^max 

NA 

45 

96.138 

91 


f 

NA 

6 

6 

6 


T 

NA 

37 

46.878 

45 



Table 4.8: Comparison of actual parameters and path-finding model results for one target mul¬ 
tiple leg scenario. 


4.7 Results 

The core effort of this thesis concerns the need to provide target parameters to the optimization 
model in Section 2.3. Thus far, we demonstrate the capability to take intelligence data concern¬ 
ing departure and arrival locations, speed, and departure time to produce a probabilistic heat 
map over time. With the heat map as input, we show the ability to determine the necessary 
target parameters to build a search vehicle allocation plan. This section briefly describes the 
results of the optimization model with our input parameters verse the results using the actual 
target parameters for several target scenarios. 
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Final Model Fit to Expected Location for 
One Target on a Multiple Leg Path 




Case 6 



Figure 4.26: Output of the path-finding model for a multiple leg scenario, showing the target’s 
track through the AOR and the target’s track in both the x and y directions with respect to time. 


4.7.1 Optimization Results 

Consider a scenario with two targets, where target 1 is an SPSS in the eastern Pacific and target 
2 is a go-fast in the Caribbean Sea. Target 1 follows a track consisting of two legs and target two 
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is on a track with a single leg. Figure 4.27 depiets the aetual track and the traek approximated 
by the path-finding model. 



Figure 4.27: Case seenario with two targets where target 1 is a SPSS operating in the eastern 
Pacifie and target 2 is a go-fast operating in the Caribbean Sea. 


Figure 4.28 depiets the optimal seareh plans based on the aetual traek and the results of our 
path-finding model. In both oases the searcher leaves from home base and travels to seareh for 
target 2 and then to seareh for target 1. The searoher’s track follows the points in sequential 
order, starting and ending at home. In eaoh instance the seareh locations and times of the 
search are nearly identioal, in addition to the expected amount of drugs found. The fact that 
the two seareh plans are so similar suggests that our path-finding model adequately determines 
the target’s track since we are comparing optimal search plans with the parameters input into 
the enhanced probability model and the output from our path-finding model. This example 
illustrates a key link between our probability and path-finding models with the seareh vehicle 
optimization model. 


58 






Figure 4.28: The optimal search plan based on the results of the path-finding model and the 
actual track. The two search plans are nearly identical. The searcher leaves from home base 
and travels to search for target 2 and the to search for target 1. The searcher’s track follows the 
points in sequential order starting and ending at home. The green track represents the search 
plan based on the true intelligence inputs, and the track in orange represents the search plan 
based on our path-finding model. 
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CHAPTER 5: 
Conclusion 


Illicit drug-trafficking is a major concern of the United States and eontinued expansion of drug- 
traffieking signifieantly affeets the national seeurity interests of the United States. This thesis 
addresses a partieular aspeet of illegal drug-traffieking in the eastern Paeifie and Caribbean Sea. 
We focus on taking intelligence information eoneerning a drug traffleker’s transit route to ereate 
a spatial temporal probability map and subsequently determine the traffieker’s most likely route 
through the region. Last, we address the aeeuraey of our results by eomparing the optimal 
seareh plan with our path-finding model as input and with the true ease parameters as input to 
the optimization model. 

We begin with a fundamental probability model to aeeount for uneertainty in the departure and 
arrival loeations of a single drug traffieker who travels along a linear path. Next we enhanee 
this model to aeeount for additional uneertainty in the traffieker’s speed and departure time. 
Additionally, we build in funetionality to handle multiple drug traffiekers, multiple eourses and 
speeds, and the ability to output results showing the expeeted amount of drugs. 

With uneertainty in the target’s transit route represented by several random variables, we ereate 
a probability map of the target’s traek over time. Using these probabilities, we apply regression 
teehniques to determine the eharaeteristies of the target’s traek needed to run a seareh vehiele 
optimization model. The optimization model requires inputs for eaeh segment of the target’s 
path where on eaeh segment the target is traveling along a eonstant heading, with a eonstant 
speed, and a eonstant AOU. Departure time and arrival time on eaeh segment identifies where 
one segment ends and the next begins. Segment ehanges oeeur based on ehanges in speed, 
heading, or AOU. Together these route eharaeteristies represent the inputs to the seareh vehiele 
optimization model. Our proeedure results in an optimal seareh plan suffieiently elose to one 
developed using the probability model inputs. 

We illustrate the models’ funetionality with two ease studies of varying eomplexity. First, a 
simple straight path between the origin and destination shows the path-finding model’s ability to 
identify the target’s starting loeation, destination, speed, and assoeiated travel times. Combining 
these route eharaeteristies, we are able to reproduee the target’s likely traek from start to finish. 
Supplying our path-finding results into the seareh vehiele alloeation model shows similar seareh 
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plan results. Next, we demonstrate a more eomplex seenario with a target traveling to the 
destination with a eourse ehange during the journey. Our model aeeurately identifies this eourse 
ehange and reproduees the target’s traek. Again, the resulting seareh plans are quite similar. 
These two ease studies are representative of aetual drug traffieker routes but additional analysis 
is needed to validate the model’s ability. 

Our probability model and path-finding model establish additional areas of researeh to aid 
JIATF South in their mission. We address the issue surrounding the AOU over time with three 
different methods where the ehosen method involves a greedy algorithm. While this method re¬ 
sults in areas mueh smaller than the other two approaehes, our overall approaeh was relatively 
informal and further researeh eould apply rigorous mathematieal methods to the AOU problem. 

We examine several eases with various target eharaeteristies and eompare our path-finding 
model results to the original probability model inputs. Our limited eomparisons demonstrate 
promising results but a more formal eomparison is neeessary to uneover any modeling shortfalls 
that may be present. Additionally, we eompare a few seareh vehiele plans using our path-finding 
model and the original ease inputs. Again, a more thorough investigation should be undertaken, 
possibly involving design of experiments or other faetors to determine where the model might 
not perform as well. 

We do not address eases with a large departure uneertainty. These are the eases referred to as 
"ehannel" oases in this thesis. They eonsist of seenarios where a target has a large departure 
window, on the order of ±24 hours or more. The probability model ereates a long ehannel 
of probability that grows from the origin and reaehes the destination, effeotively oonneeting 
the two looations. In other words, the probability mass is never eompletely on the water and 
the target oan be anywhere along the ehannel from the earliest departure time to the latest 
arrival time. The optimization model handles these eases in a different manner and requires 
different inputs, thus follow-on analysis to develop different algorithms eonsidering these eases 
is neeessary. 

Our probability model does not aeeount for the impaet of weather on the drug traffieker’s move¬ 
ments. A eomparison of the path-finding model results utilizing input from our model and one 
with weather eharaeteristies would highlight whether or not those items are important faetors 
for eonsideration. 
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